Spinor Bose-Einstein condensate flow past an obstacle 
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We study the flow of a spinor (F = 1) Bose-Einstein condensate in the presence of an obstacle. We 
consider the cases of ferromagnetic and polar spin-dependent interactions and find that the system 
demonstrates two speeds of sound that are identified analytically. Numerical simulations reveal the 
nucleation of macroscopic nonlinear structures, such as dark solitons and vortex-antivortex pairs, 
as well as vortex rings in one- and higher-dimensional settings respectively, when a localized defect 
(e.g., a blue-detuned laser beam) is dragged through the spinor condensate at a speed larger than 
the second critical speed. 
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I. INTRODUCTION 

Over the last decade, we have seen an enormous growth 
of interest and a related diversification of the physics of 
atomic Bose-Einstein condensates (BECs) Q . A sig- 
nificant aspect of this ever-expanding interest is the in- 
tense study of macroscopic nonlinear excitations, such as 
solitons and vortices, which can arise in BECs 0. In 
fact, the emergence of such macroscopic coherent struc- 
tures in the many-body state of the system establishes 
a close connection between BECs and other branches of 
physics, such as e.g., optics and the physics of nonlin- 
ear waves. Within this interface of atomic and nonlinear 
wave physics, in recent years there has been an increas- 
ing focus on the study of multi-component BECs [3], 
and particularly spinor condensates [3, l5|. The latter 
have been realized with the help of far-off-resonant op- 
tical techniques for trapping ultracold atomic gases Q 
which, in turn, allowed the spin degree of freedom to be 
explored (previously frozen in magnetic traps). This rel- 
atively recent development has given rise to a wealth of 
multi-component phenomena, including the formation of 
spin domains [7] and spin textures 1^, spin- mixing dy- 
namics [9.], dynamic fra gme ntation [1^, and the dynam- 
ics of quantum phases At the same time, macro- 
scopic nonlinear structures that may arise in spinor BECs 
have also been investigated. Such structures include 
bright [13, [3, El, dark [H, and gap solitons [3, as 
well as more elaborate complexes, such as bright-dark 
solitons and domain- walls [18]. 

A relevant direction that has been of particular inter- 
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est concerns the study of the breakdown of superfiuidity 
and the concomitant generation of excitations in BECs. 
In particular, much experimental and theoretical effort 
has been devoted to the understanding of important rel- 
evant concepts, such as the critical velocity introduced 
by Landau, sound waves and the speed of sound, and the 
emergence of vortices and solitons [i, ^2] . From the the- 
oretical point of view, the Gross-Pitaevskii (GP) equa- 
tion has been used to study the fiow of a BEC around 
an obstacle or, equivalently, the effect of dragging a lo- 
calized potential (such as a blue-detuned laser beam) 
through a BEC. In this context, it has been predicted 
theoretically IIO, Hil, |22] and recently observed ex- 
perimentally [23] that when the speed of the "localized 
defect" exceeds a critical speed, then dark solitons are 
formed in quasi one-dimensional (ID) condensates. On 
the other hand, in higher-dimensional [e.g., quasi two- 
dimensional (2D)] settings, theoretical studies [24] have 
shown that a similar procedure leads to the formation 
of vortices (or more precisely to vortex-antivortex pairs, 
due to the conservation of total topological charge). Im- 
portantly, experimental consequences of this procedure, 
such as an onset of heating and dissipation were mon- 
itored experimentally [25[. Other relevant theoretical 
works include studies of the breakdown of superfiuidity, 
the onset of dissipation, and the associated Landau cri- 
terion [26]. More recently, dragging of an obstacle in a 
two-component BEC was studied in Ref. [27]. In this lat- 
ter study it was established that two distinct "speeds of 
sound" arise and the form of the ensuing nonlinear struc- 
tures (e.g., dark-dark or dark-anti-dark soliton pairs in 
ID, and vortex- vortex or vortex-lump pairs in 2D) de- 
pend on how the value of the obstacle speed compares to 
the values of the critical speeds. 

In this paper, we consider the dragging of a localized 
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defect through an F = 1 spinor condensate with repul- 
sive spin-independent interactions and either ferromag- 
netic or anti-ferromagnetic (polar) spin-dependent inter- 
actions. In the framework of mean-field theory, this sys- 
tem is described by a set of coupled GP equations for 
the wavefunctions of the three hyperfine components. A 
key question that emerges in this F = 1 spinor BEG set- 
ting is how many critical speeds may be available. A 
naive count based on the three-component nature of the 
system (and by analogy to the two-component setting 
bearing two such critical speeds) would suggest the pos- 
sibility for three distinct critical speeds. However, as we 
illustrate below, an explicit calculation reveals that there 
exist only two such critical speeds in the system, due 
to the particular nature of the nonlinear ity. Moreover, 
our numerical simulations illustrate that the crossing of 
the lower of the two critical speeds does not appear to 
lead to the formation of nonlinear excitations. On the 
other hand, for defect speeds larger than the second crit- 
ical speed, our simulations illustrate that dark solitons 
emerge in the ID setting, vortex-antivortex pairs in the 
2D setting, and vortex rings are shown to arise in the 
fully three-dimensional (3D) setting. 

The presentation of our results is structured as fol- 
lows. In Sec. II, we develop an analytical approach for 
computing the relevant critical speeds, by generalizing to 
the spinor setting the arguments of Ref. [19]. Then, in 
Sec. Ill, we numerically test the relevant predictions in 
ID, 2D and 3D settings. Finally, in Sec. IV, we summa- 
rize our findings and point to some important remaining 
questions along this vein of research. 

II. MODEL AND ITS ANALYSIS 

In our analytical approach, we will consider a quasi- 
ID spinor F = 1 BEG with repulsive spin-independent 
interactions. In the framework of mean-field theory, this 
system can be described by the following normalized GP 
equations flTl. fTsj: 

idt^±i = i^oV^±i+r[(|V^±ip + |V^o|'-|V^Til')V^±i] 

+rv^oV^;i, (1) 

+2rV^_iV^o>+i, (2) 

where Hq = -(1/2)S^ + V{x',t) + ntot, while ntot = 
+ IV^oP + 1^+1 P is the total density and V{x;t) 
is the external potential. The latter is assumed to take 
the following form: 

V{x; t) = ^^^^^ + Vo exp{-a{x - stf). (3) 

The first term in the right-hand side of Eq. (|3]) repre- 
sents a typical harmonic trapping potential of normal- 
ized strength while the second term accounts for a 
localized repulsive potential (e.g., a blue-detuned laser 



beam), of strength Vq and width a~-^, that is dragged 
through the condensate at speed s. Note that our an- 
alytical results will be obtained below for the case of 

= (which still contains the fundamental phenomenol- 
ogy), but were also tested in the numerical simulations 
for 1] / (and were found to persist in the latter case). 
Finally, the parameter r in Eqs. ([I])-© expresses the 
normalized spin-dependent interaction strength defined 
as r = (a2 — ao)/(ao + 2a2), where ao and a2 are the 5- 
wave scattering lengths in the symmetric channels with 
total spin of the colliding atoms F = and F = 2, re- 
spectively. Note that r < and r > correspond, re- 
spectively, to ferromagnetic and polar spinor BEGs. In 
the relevant cases of ^^Rb and ^^Na atoms with F = 1, 
this parameter takes values r = —4.66 x 10~^ [28] and 
r = +3.14 X 10~^ [29], respectively, i.e., in either case, it 
is a small parameter in Eqs. ([I])-©. 

We now seek uniform stationary solutions of the GP 
Eqs. (HI)-® (with V = 0) in the form 

tjji = Aexp(— z/i+it) exp(z^+i), 

= Cexp(— z/i_it) exp(z^_i), 

where A, 5, C and Oj (with j G { — 1,0,+!}) represent, 
respectively, the amplitudes and phases of the hyperfine 
components, and /nj are their chemical potentials. In our 
analysis below we will assume that A 0^ B 0^ C 0, 
as that will provide us with genuinely spinor (i.e., three- 
component) states; otherwise the system is reduced to a 
lower number of components. In fact, the analysis for the 
one-component case has been carried out in Ref. [l^, 
while in the two-component case considerations analo- 
gous to the ones that we will present below have been 
put forth in Ref. Under the above genuinely three- 

component assumption, we substitute the stationary so- 
lutions into the GP Eqs. ([I])-© and obtain the following 
set of equations: 

B'^C 

/i+i = ntot + r(A2+52-C^)+pr— ^, 
/io = ntot + r(A2+C^)+2prAC, 

B'^A 

/i-i = ntot^r{C^ ^B^ -A^)^pr-^, 

where ntot = + 5^ + C^. 

In the above expressions phase matching conditions 
were used, as is usual when one has parametric inter- 
actions: these read 2/io = M+i + M-i for the chemical 
potentials, and AO = 26o — (^+i + 6>_i) = or tt for the 
relative phase between the hyperfine components [l^, [sO] • 
The factor p = ±1 on the last term of each of the above 
equations results from considering = or tt, respec- 
tively. In the case where the three chemical potentials 
jiij are different, it can be found that it is not possible to 
satisfy the above assumption that each of the amplitudes 
A, 5, C should be nonzero. Hence, we will hereafter 
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focus on the case of /i+i = /io = M-i = In the lat- 
ter case, it is straightforward to algebraicahy manipulate 
the equations and find that there exist only two classes of 
possible stationary solutions with a free parameter (for a 
given /i). These solutions are as follows: 

A = -pC, B = ±V/i-2C2, /i > 2C^ (4) 
A = -pC ± ]Jy^^ ^ = /i > 1 + r. (5) 

and are simply the "anti-phase-matched" and "phase- 



matched" type solutions of Ref. respectively. Note 
that in addition to these solutions, there exists another 
one with no free parameters (i.e., all amplitudes are di- 
rectly dependent on /i) which is a particular case of Eq. 
([5]), namely 



A = pC, B = ±V2C^, /i = 4(l + r)C^ (6) 

Let us now consider the GP Eqs. ([I])-© looking for 
stationary solutions in that frame the GP equations be- 
come: 



-isda:tlj±l = -^dlllj±i +ntot^±l +r(|7/;±ip + iV^oP - |V^TlP)V^±l + ^V^O^Tl ~ '^^il' 



+1 



(7) 
(8) 



where we have slightly abused the notation by replacing 
the traveling wave variable £^ = x — st with x for simplic- 
ity. We now decompose the amplitudes Rj and phases (j)j 
of the order parameters according to ipj = Rj ex.p{i(j)j)^ 
and impose the phase matching condition ^i+^-i = 2^o, 
to obtain the following equations: 



(9) 



d^h = s(i- i—y 



dlRi = -s^ Ri 



Rl 



dlRo = -s^Ro 



2ntoti?i (10) 
2r{Rl + i?o - R-i)Ri + 2prRlR-i - 2/ui?i, 

(11) 
(12) 



a,0o = s{i- i—f 



2 / I 2ntot^o 



dx(l)-i 

dlR-i = -s^ [ R-i 



2r{Rl + Rii)Ro + 4.prRoRiR-i - 2^Rq, 



(13) 



2ntot^-i 



(14) 



+ 2r{R\ -^Rl- R\)R-i + 2prRlRi - 2/ii?_i. 

Notice that in these equations, the asymptotic states A, 
5, and C arise naturally due to the integration of the 
equations for the phases. 

In order to seek instabilities of the steady state fiow 
at different fluid speeds, we now linearize around the 
asymptotic states, according to Ri = A + eri(x), Rq = 
B + er^ix) and R-i = C + er_i(x) (where e is a formal 
small parameter). Substituting the above expressions 



r 



into Eqs. ([9|)- (p!5|) . we obtain a system of three second- 
order ordinary differential equations; the latter, can be 
readily expressed as a system of six flrst-order equations 
of the following form: 



d_ 

dx 



where r'^ = drj/dx and 
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M= {rriij} 



and the non-zero matrix elements of M are given by the 
following expressions: 

m2i = -45^ + 2(3^2 + 52 + 0^) + 
2r{3A^ ^ B^ - C^) - 2/1, 

77723 rn^i 4:AB{1 + r) + AprBC, 

m25 m6i=4AC{l-r)^2prB'^, 

m43 = -452 + 2(352^^2 + 02) + 
2r(A2 + c2)+4prAC-2/i, 

77145 rriQs 4BC{1 + r) + 4prAB, 

^65 = -452 + 2(3^2+52+^2) + 

2r(3C2 + 52-A2)-2/i. 
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Notice that in order to derive the above system of ordi- 
nary differential equations, we have partially simplified 
the problem, assuming no perturbations in the phases. 
In such a more general case, however, the full first-order 
ODE system incorporating phase perturbations is in fact 
twelve-dimensional and is not analytically tractable. We 
have found (not treated explicitly here) that this system 
can be analyzed only in some special cases, such as 5 = 0, 
yielding the same results for the critical defect speeds, as 
will be presented below (see also the discussion of Section 
II.B in Ref. 

We now follow the approach used for the one- 
component GP equation in Ref. [l^ (subsequently gen- 
eralized in the two-component case in Ref. [23]) to deter- 
mine the critical speeds of the defect. In particular, the 
critical speeds can readily be found upon computing the 
eigenvalues of the matrix M and setting them equal to 
zero. The violation of this threshold condition is tanta- 
mount to the emergence of a number of oscillatory modes 
that enforce too many constraints and prevent the exis- 
tence of localized solutions for a generic obstacle poten- 
tial, as is explained in detail in Ref. [l^. It is straightfor- 
ward to examine this condition both in the case of p = 1 
and of p = —1. We will demonstrate below the case of 
p = 1 for definiteness. In this case, by considering the 
stationary state of the form of Eq. we obtain two 
different speeds of sound, namely, 

ci ^/rc2, C2 = v7^. (15) 

It is clear that the first critical velocity is characteristic 
for the spinor F = 1 condensate under consideration (as 
it depends on the normalized spin-dependent interaction 
strength r), while the second one is the standard speed of 
sound appearing in the one-component GP equation [19] 
(note that this speed may also appear in the case of the 
spinor F = 1 BEG as well pJTj). It is worthwhile to point 
out that as r > for anti-ferromagnetic bosonic spin- 
1 atoms (e.g., ^^Na), while it is r < for ferromagnetic 
ones (e.g., ^^Rb), then the first speed of sound is relevant 
(i.e., will only exist) in the case of, e.g., the polar ^^Na 
spinor condensate. On the other hand, by selecting the 
stationary states of the form of Eq. ([5]), then again we 
find two critical speeds, which are now given by: 



ci = Y YTr^'^' " ^^^^ 

In this case, it is clear that the first critical speed will ex- 
ist only in the ferromagnetic spinor BEGs (such as ^^Rb), 
but not in anti- ferromagnetic ones (such as ^^Na); nev- 
ertheless, it should be noted that since the normalized 
spin-dependent interaction strength is small in both cases 
of ^^Rb and ^^Na condensates (r = O(10~^) as discussed 
above), the lower critical speeds are approximately the 
same. 

We now test these analytical predictions by drag- 
ging a localized defect (e.g., a blue-detuned laser beam) 
through the condensate at different speeds characteriz- 



ing the three regimes, namely (a) < 5 < ci < C2, (b) 
Ci < 5 < C2, and (c) Ci < C2 < s. 
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FIG. 1: Dependence of the critical velocity on the defect 
strength Vq. The result shown corresponds to the station- 
ary solution of Eq. {g]) with A = -0.5 = -C and B = 
(/i = 1). It was confirmed that, e.g., for the solutions of 
Eq. (|5]), the analytical and numerical results were indistin- 
guishable up to three decimal places. 



III. RESULTS OF NUMERICAL STUDIES 

A. One-dimensional setting 

Our more detailed results concern the ID setting, 
where we explore the full two-parameter space of speeds 
s and defect strengths Vb, for a = 2 in the case of 
the anti- ferromagnetic ^^Na spinor BEG, characterized 
by the spin-dependent interaction strength r = 0.0314. 
Figure [1] illustrates the threshold above which coherent 
localized excitations are emitted from the defect as it 
propagates through the condensate. A typical example 
of the evolution process for (a) < 5 = 0.025 < ci < C2, 
(b) ci < s = 0.325 < 02 and (c) ci < C2 < s = 0.335 is 
shown in Fig. [2l Several comments are in order here: 

• As expected from the analytical predictions, when 
the defect speed is below both critical values ci = 
^/JIr^ and C2 = ^/fl^ the defect moves through the 
atomic cloud without emission of any nonlinear ex- 
citation. An oscillatory structure is radiated at the 
initial time (similarly to what has been observed 
earlier, e.g., in Ref. [27]) both at the front, as well 
as at the rear of the defect, moving with the speed 
of sound; however, no further such radiation is ob- 
served. 

• Remarkably, for speeds intermediate between c\ 
and C2, we do not observe any modification in 
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FIG. 2: (Color online) Time evolution of the three hyperfine 
components (top, middle and bottom rows, respectively) for 
different defect velocities: s = 0.025 < ci (left column), a < 
s — 0.325 < C2 (middle column) and s — 0.335 > C2 (right 
column); here, ci = 0.058 and C2 = 0.329 are the two critical 
velocities, while the defect strength is taken to be Vo = 0.9. 
The initial condition corresponds to the stationary state with 
wavefunction amplitudes given in Eq. (|4]). While the results in 
the left and middle columns show a steady flow (apart from an 
oscillatory structure that is detached at t = 0) the evolution 
shown in the panels of the right column is characterized by 
the emission of dark solitons, even from the early stages of the 
process. Notice that the analytically predicted first (lower) 
critical velocity would fall between the velocities of the results 
depicted in the left and middle columns but no significant 
change is observed in the dynamics between these two cases. 



the dynamics. This means that the first criti- 
cal speed ci does not appear to be activated by 
the system. This finding is even more surprising 
in light of the fact that for r < 0, this critical 
speed has been recognized to be directly connected 
to the quasi-momentum (wavenumber) associated 
with the modulational instability of the ferromag- 
netic spinor condensate [30] (see also the relevant 
discussion in Ref. [17.]). Nevertheless, in all of our 
simulations, both in ID and in higher dimensions, 
we have definitively confirmed the apparent phys- 
ical irrelevance of this first critical speed (which 
is the lower nontrivial critical speed in the spinor 
BEC case). It should be noted here that this same 
feature has been confirmed for cases where the spin- 
dependent interaction strength r was artificially in- 
creased to considerably larger values (by an order 
of magnitude in comparison with its physically rel- 
evant value of r = 3.14 X 10"^ for ^^Na). 

• When the defect speed is larger than both criti- 
cal ones, there is a clear emission of dark (in fact, 
gray) solitons, which travel in a direction opposite 
to that defect, with velocities less than the speed of 
sound. Similarly to one- [H, [lO, Hll, [22[ and two- 
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FIG. 3: (Color online) Vortex and ant i- vortex pairs nucleated 
by a moving defect in the 2D spinor condensate. The left col- 
umn shows the density of the spinor condensate's component 
-00 at different times (indicated in the panels). The location 
and extent of the moving defect is depicted by the oval line 
corresponding to an iso-contour of its strength at 10% of its 
maximum. The right column shows the corresponding vortic- 
ity (defined in the text), clearly illustrating the presence of 
vortex and anti-vortex pairs. The case depicted here corre- 
sponds to k;^ = 8, a = 2, and Vb — 0.9 while the defect speed 
is taken to be s = 0.6 > C2. 



component [27[ settings, the solitons temporarily 
"alleviate" the super-critical nature of the flow, but 
eventually they are separated enough from the de- 
fect that another such excitation emerges. For this 
reason, the emission seems to be regularly spaced 
as shown in the right panels of Fig. [2l 

• As the strength of the defect ^ 0, the critical 
speed C2 observed from the numerical simulations 
tends asymptotically to the one theoretically pre- 
dicted from the analysis above, i.e., ^JJl = 1 in this 
case (this feature has been confirmed for different 
values of /i, such as /i = 2 and /i = 4). However, 
similarly to what was observed in Refs. [l9|, [27^ . 
as the strength of the defect increases, the value 
of the critical speed accordingly decreases (since 
nucleation of dark solitons is easier for the lower 
density BEC). 

• Finally, we note that in all our simulations (even 
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FIG. 4: (Color online) Evolution of the maximum (dark/red) 
and minimum (light /yellow) of the vorticity iso-contours in 
(x,y,t). One can clearly discern the emergence of vortex (in 
dark/red) and antivortex (in light /yellow) pairs, as the defect 
moves along the x- direct ion. The left and right panels corre- 
spond, respectively, to the cases depicted in Figs. [3] {wy — 8) 
and[5](^^ 20). 



in higher dimensions, see below) the three spinor 
components were locked to each other through 
IV^ilVA^ ^ \M^/B^ ^ IV^-ilVC^. This tight 
restriction is presumably related to the fact that 
we only observed one critical nucleation speed in 
our simulations: as all the components are tightly 
locked to each other, they behave like a single com- 
ponent and thus only one critical speed is observed. 



B. Two-dimensional setting 

In the 2D case, motivated by the recent work of Ref. 
ps] in the single component case, we consider a defect 
which is localized along the x-axis (with a width a~^) 
but elongated along the ?/-axis (with a width Wy > a~^), 
namely: 



V 



— exp [—a[x 



tanh 



stf 



1 



tanh 



1 



Once again, we find that (i) no emission of nonlinear ex- 
citations is present for s < C2 and that (ii) the emission of 
nonlinear excitations, which now have the form of vortex- 
ant ivortex pairs, arises for speeds larger than the critical 
speed C2. Figure [3] illustrates the case of s = 0.6 > C2 
for the defect strength Vq = 0.9, while the defect width 
along the ^/-direction is il'^ = 8. In addition to show- 
ing the density, the figure shows the vorticity defined as 
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FIG. 5: (Color online) Same as Fig. [3] but for a wider de- 
fect with Wy — 2^. One can clearly observe the formation of 
numerous vortex anti- vortex pairs in the wake of the defect. 



V X 'Uf, where the fluid velocity is given by 



(17) 



for a given hyperfine-component ip. In the contour plots 
of the vorticity co, the emergence of vortex- antivortex 
pairs is immediately evident in the supercritical case 
shown. In fact, in order to provide a more clear sense of 
the temporal dynamics and the nucleation of the coher- 
ent structure pairs, we show in the left panel of Fig.[4]the 
spatiotemporal evolution of the iso-contours of vorticity 
for the same numerical simulation as in Fig. [3l The vor- 
ticity renders transparent the emergence of the different 
vortex pairs at different moments in time (and accord- 
ingly different locations in x, as the defect travels). 
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Finally, in Fig. [5l we also show a case example of a 
considerably wider defect, with a width Wy = 20. It can 
be seen that in this setting, the region of low density 
caused by the defect is far wider, in turn leading to a 
breakup into a large number of vortex pairs that can be 
identified not only by the density minima, but also even 
more clearly (including their topological charge) by the 
vorticity panels. In the right panel of Fig. [H we depict 
the corresponding spatiotemporal evolution of the vor- 
ticity. We note that, as it was the case in the ID setting, 
all our simulations suggest that the three spinor com- 
ponents remain essentially locked satisfying the relation 
IV^ilV^^ ^ IV^olV^^ ^ For this reason, we 

only depict the dynamics of the tpo component in all of 
our results. 



C. Three-dimensional setting 

Finally, we also performed 3D simulations, using a 
three-dimensional generalization of the potential, elon- 
gated along the z-direction, namely 

^0 



V 



16 

tanh 



exp (— a(x 

( 



sty 



y 



^1 + 1 



tanh 



1 



tanh 



1 



1 



X tanh (^z 

with a = 2, Wy = 8, Wz = 4 and Vq = 0.9. In Fig. [6] 
we depict the results for s = 0.8 (i.e., above the sec- 
ond critical speed). Panels (a)-(e) depict the iso-density 
contours for |?/^oP at different times while panel (f) de- 
picts a superposition of the isocontours for the norm of 
the vorticity field. As it can be seen from the figure, 
a vortex ring is formed in the 3D spinor condensate, as 
a result of the supercritical nature of the chosen speed 
s = 0.8 > C2. The iso-density contours of IV^oP clearly 
show a depletion of atoms around the vortex ring that 
is nucleated in the wake of the defect-induced region of 
density minima. It is worth stressing again that the dy- 
namics of the different components seems to be locked 
such that l^ipM^ ^ l^olV^^ ^ \i^-i?/C^ and, there- 
fore, we only depict results for ipo in the 3D case as well. 

We also performed simulations for larger defect speeds 
giving rise to a rich and complex scenario of multiple vor- 
tex rings nucleations, collisions, collapses and splitting. 
A typical case is shown in Fig. [3 that corresponds to the 
same parameters as in Fig.[6]but for a larger defect speed 
{s = 1). The main characteristics of the evolution can 
be summarized as follows: a first vortex ring is nucle- 
ated in the wake of the defect at about t = 12; around 
t = 17 a second vortex ring is nucleated while the first 
vortex ring starts to shrink until it eventually disappears 
around t = 23; shortly after this, the second ring deforms 
and splits into two separate vortex rings around t = 30. 
IV. CONCLUSIONS AND FUTURE 
CHALLENGES 

We have studied the motion of a localized defect 
through a spinor F = 1 condensate. Despite the three- 



component nature of the system, our systematic analysis 
of the small-amplitude excitation problem revealed that 
the nature of the nonlinearity is such that there appear 
not three, but merely two critical speeds in the system; 
these were identified analytically for the families of sta- 
tionary uniform states of the system. Our numerical sim- 
ulations tested the dynamics for different values of defect 
speeds in comparison to the two critical speeds. Surpris- 
ingly, it was found that the lower one among the two 
critical speeds is not activated and no emission of non- 
linear wave excitations emerges when this threshold is 
crossed. On the other hand, when the defect speed ex- 
ceeds the second critical one, then emission of coherent 
structures arises independently of dimension; the result- 
ing waveforms are dark solitons in the one-dimensional 
setting, vortex-antivortex pairs in two dimensions, and 
spinor vortex rings in the fully three-dimensional case. 

While the present study showcases an experimentally 
accessible mechanism for producing nonlinear excitations 
in spinor BECs, a number of interesting questions are still 
outlying. In particular, perhaps the most relevant ques- 
tion from a theoretical point of view involves acquiring 
a full understanding of why the first critical speed does 
not seem to be explored by the system, contrary to what 
might be expected from the two-component case ana- 
lyzed in Ref. [27]. Another direction of potential interest 
could be to explore in this multi-component system what 
would happen if the speed of the defect becomes consid- 
erably larger than the critical speed, in which case, and 
in the one-component setting, a convective stabilization 
of oblique dark solitons has been reported Also in- 
teresting would be to study the formation of shock waves 
[33] in the spinor systems. Work along these directions 
is currently in progress and will be reported in future 
publications. 
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FIG. 6: (Color online) Vortex ring formation in the supercritical 3D setting. Panels (a)-(e) depict iso-density contours 
corresponding to |'0o|^ = 0.3 at the indicated times. Note that the extent of the moving defect is clearly visible in these panels 
(it corresponds to the rightmost flat oval shape that is created by the atomic density depletion due to its presence). Panel (f) 
shows typical isocontours of the norm of the vorticity of ipo at times t = 12,14, ... , 30 (left to right). In this case we use a 
defect with Wy = 8, Wz = 4, a = 2 and speed s — 0.8 > C2. 
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FIG. 7: (Color online) Similar to Fig. [6] but for a slightly larger defect speed s — 1 > c^. The top row of panels depicts the 
iso-density contours of V^o while the bottom row depicts the respective vorticity isocontours. Note the successive nucleation of 
two vortex rings. The first one shrinks and collapses into itself between t — 22 and t — 23, while the second ring deforms and 
eventually splits into two separate vortex rings between t = 29 and t = 30. 
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